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Abstract 

We discuss the end-point behavior of the pion distribution amplitude (DA) and calculate 
its slope using QCD sum rules with nonlocal condensates. This is done in terms of the 
standard derivative and also with the help of an "integral derivative" , recently obtained by 
us. Our approach favors a value of the slope of the order (or less) of the asymptotic DA 
and is in clear disagreement with flat-type pion DAs. 

1 Introduction 

The main technique to analyze hard exclusive processes within QCD, is provided by the factor- 
ization of the underlying dynamics into a hard and a soft part. The hard part forms the partonic 
amplitude of the subprocess at a large value of the momentum transfer and is amenable to QCD 
perturbation theory. The soft part depends on the distribution amplitude of the hadron(s) and 
contains the dynamics at typical hadronic scales; it has, therefore, to be determined by nonper- 
turbative methods or be extracted from experimental data. The collinear factorization applied 
to the transition form factor (FF) of two far off-shell photons to a pion leads to the convolution 
of these two parts, which, at the leading order of twist two, reads (x = 1 — x) 



F^{Q\ g 2 ) = ^-U £ Q2 - \ q2x M*) dx + 0(l/Q 4 



(1) 



modulo twist-four terms, ignored here. The main ingredient of the above equation is the pion 
DA (Ptt(x) which encodes all unknown binding effects of the pion state. At the considered level 
of twist two, it is defined by the following universal matrix element [TJ 



(0\d(zh^ 5 [z,0]u(0)\n(P))\ z 2 =0 = iUP» f e ix ^ipt 2) (x,^)dx 

Jo 



(2) 



where x is the longitudinal momentum fraction carried by the valence quark (x for the antiquark) 
in the pion and the path-ordered exponential, i.e., the light-like gauge link, 



[z, 0]=V exp 



ig / t a Al{y)dyK 



o 



(3) 
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ensures gauge invariance. It is useful to expand the pion DA in terms of the Gegenbauer har- 
monies 6xx C n (2x — 1) which provide the eigenfunctions of the one-loop Efremov-Radyushkin- 
Brodsky-Lepage (ERBL for short) evolution equation [21(3]. One finds for c^( t=2 ) at the typical 
hadronic scale 



^ t=2 \x- /xg) = <p***{x) \l + a 2 (fxl) Cl /2 {2x - 1) + a 4 ( / Uo) C 3 J 2 (2x - 1) 



(4) 



in which the asymptotic (abbreviated by Asy) pion DA appears: 9? Asy (x) = 6xx. By virtue 
of the leptonic decay it — > /U + z^, one obtains the normalization tp$~ 2 '(x, /U 2 ,) dx = 1, which 
fixes ao = 1. 

While a process involving two photons with large virtualities is theoretically preferable, be- 
cause one can safely apply QCD perturbation theory, experimentally, the asymmetric kinematic 
with one of the photons being quasi real is more accessible. Indeed, such measurements have 
been carried out by several collaborations, namely, the CELLO 0], the CLEO [5], and, most 
recently, the BaBar Collaboration [6j. Taking the limit q 2 — > in convolution ([1]), one finds 
that the FF for the 7*7 — > 7r° transition is actually given by the inverse moment of the pion 
DA 



(^7r(x) dx . 



(5) 



Therefore, this quantity is one of the key elements of the pion-photon transition FF. Because 
this form factor has such a simple structure within QCD, it has attracted over the years the 
attention of many theorists (see, e.g., [7HTT] and references cited therein). 

However, the most recent measurement of this observable by the BaBar Collaboration [6] 
has provided controversial results, because, unexpectedly, the high-energy data points above 
10 GeV 2 grow with Q 2 — see Fig. [TJ At moderate values of the momentum transfer, up to 
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Figure 1: Predictions for the photon-to-pion transition FF calculated in next-to-leading order (NLO) light- 
cone sum rules using a Breit-Wigner model for the meson resonances. The pion is parameterized in terms of the 
following DAs: Asy — dashed-doted line, BMS "bunch" [10] — shaded (green) strip, CZ model [12] — upper 
long-dashed (red) line, flat-top model [T3] — short-dashed (blue) line. Experimental data are shown on both 
panels using the following notations: BaBar data [6] — (red) boxes with error bars, CELLO data [4] — (black) 
diamonds with error bars, CLEO data 5 — (violet) triangles with error bars. The horizontal solid line denotes 
the asymptotic QCD prediction %/2/vr. The left panel provides a zoom-in view into the low and intermediate 
range of Q 2 , whereas the right panel shows the results for the whole interval of momenta probed by the BaBar 
experiment. 

9 GeV 2 (see the left panel of Fig. [1]), the new high-precision BaBar data agree well with the 
previous CLEO data [5]. From the second column of Table [T] and the left panel of Fig. [H we 
may conclude that all data up to 9 GeV 2 can be best described by pion DAs that have their 
end-points strongly suppressed [TTJE]. A characteristic example of such a DA is provided 
by the Bakulev-Mikhailov-Stefanis (BMS) model [10] , which has been derived from QCD sum 
rules (SR)s with nonlocal condensates (NLC)s, originally developed in [15]. In contrast, the 
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high-C} 2 BaBar data (see the right panel of Fig. [T]) show, as already mentioned, an unexpected 
growth with Q 2 which cannot be understood on the basis of the collinear factorization and 
calls (see the fifth column of Table [1])) for pion DAs that have instead their end-points strongly 
enhanced |131I16|. This intriguing behavior has triggered the use of flat-type pion DAs and 
provided the main motivation for our analysis in [IT], on which we report here. The result 
obtained by Radyushkin [TO] with the flat-top model (p n (x) = 1) is shown in Fig. [T]as a dashed 
(blue) line in comparison with the BMS bunch (green strip) |10| and the CZ pion DA model 
(long-dashed line in red color) [IS]. As we see from this figure, the high-C} 2 BaBar data are 
rather well-described in Radyuskin's approach [13]. A recent independent analysis [18] comes 
to the conclusion that a flat pion DA, when used in a fully consistent way, yields to predictions 
for the pion's electromagnetic and transition form factors which are in striking disagreement 
with experiment. 

Table 1: Deviations of theoretical predictions for the quantity Q 2 F 1 J7r (Q 2 ) in terms of 
(ndf = number of degrees of freedom) for several pion model DAs: Asy, BMS model, and CZ model (more 
details in [11]). Predictions based on the flat-top model DA (I18|l . discussed in |13] . are also included. The second 
column shows the results for the combined sets of the CLEO [5] and the CELLO [I] data above 1 GeV 2 . The 
third column compares model predictions with all data in the interval [1, 9] GeV 2 , while the fourth column takes 
into account only the BaBar data above 9 GeV 2 . The last two columns show the values of the inverse moment 
and the standard derivative of selected pion DAs at the origin. 



DA Model 


CLEO&CELLO 


All data Q2<9 GeV 2 


BaBar g2>9 GeV 2 




¥4(0) 


Asy 


2.7 


6.35 


18.9 


3 


6 


BMS [TO] 


0.56 


0.86 


11.7 


3.15 


1.12 


cz p2] 


5.9 


33.9 


7.9 


4.5 


26 


DA [TO] 


4.15 


4.15 


1.0 







In our recent investigation [17], we revisited the QCD SR approach of [TO] focusing our 
attention on the behavior of the pion DA in the end-point region x ~ with the aim to 
understand the fine structure of the pion DA in this region vs the ansatz for the quark-virtuality 
distribution in the nonperturbative QCD vacuum. As we shall explain in Sec. [31 QCD SRs were 
mainly developed with the purpose to study the integral characteristics of the pion DA. To 
overcome this restriction, we constructed in [T7] an operator for integral derivatives of the pion 
DA. The results obtained this way supplement those found with SRs which employ the standard 
derivative of the pion DA. 

Our presentation will concentrate on the following issues, organized in sections. In Sec. [21 
we discuss the QCD SR approach with NLCs and focus on the peculiarities of the pion DA in 
the end-point region x ~ 0. In Sec. we define and study the integral derivative of the pion 
DA, the purpose being to overcome the restrictions inherent in QCD SRs which were created in 
order to probe the integral characteristics of the pion DA. In Sec. E] we also study the SRs for 
the standard derivative of the pion DA at the origin and complete our presentation by drawing 
our conclusions in Sec. [SJ 

2 Slope of the pion DA and the nonperturbative QCD vacuum 

As already mentioned in the Introduction, the fine details of the pion DA in the region around 
the origin x ~ are of crucial importance. To illustrate the differences in the end-point behavior 
of existing model DAs, we show the corresponding profiles in Fig. [2] (left panel). We depict four 
models: BMS [TO] model — solid line; CZ [12] — dashed-dotted (blue) line; flat-top DA given 
by Eq. (|18p — short-dashed (red) line; dotted line — asymptotic DA. Using the values of the 
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derivative of the pion DA, collected in the last column of Table [T] we can classify the pion DA 
models according to their end-point behavior: end-point suppressed (Asymptotic and BMS) and 
end-point enhanced (CZ and flat-top (]18p ). A zoomed-in view of the end-point characteristics 
of the above DAs is displayed in the right panel of Fig. [2] 




0.0 0.2 0.4 0.6 0.8 1.0 0.00 0.02 0.04 0.06 0.08 0.10 



Figure 2: Comparison of selected pion DA models. The left panel shows their shapes, while the right panel 
shows a magnified view of the endpoint region x ~ 0. Solid line — central line of the BMS bunch [TO]; dashed- 
dotted (blue) line — CZ model [15] : short-dashed (red) line — flat-top DA from Eq. (TT8)) with a = 0.1; dotted 
line — asymptotic DA. All DAs are normalized at the same scale fj-o ~ 1 GeV 2 . 

Relying upon the QCD SR approach with NLCs, developed in [10], we direct our attention 
to the end-point behavior of the pion DA. This will allow us to obtain predictions for the slope 
of the pion DA in the end-point region with smaller errors than those obtained from the BMS 
bunch [TU] in this region [17] . Recall that the basic idea underlying the NLC parametrization 
of the QCD vacuum, is that it has a domain structure, parameterized in terms of condensates, 
the latter possessing a certain correlation length by virtue of which the vacuum quarks ac- 
quire a non-zero average virtuality (kg) (see, for instance, [E]). To analyze the nonlocality of 
the vacuum condensate, it is useful to parameterize the lowest scalar condensate in terms o£] 
(q(0)[0, z]q(z)) = M s (z 2 ). This quantity can be related to the vacuum distribution function 
fs(a) via [E] 

oo 

M s (z 2 ) = (qq) ff s (a)e az2 ' A da, (6) 
o 

where a is the vacuum-quark virtuality. Assuming that the vacuum quarks have a fixed virtu- 
ality X 2 one has 

fs(a) = 5(a-X 2 q /2), (7) 

leading for the scalar-quark condensate [15] to the Gaussian model 

M s (z 2 ) = (q(0)q(0))e-^ s . (8) 

The parameter A 2 here represents the typical quark momentum in the QCD vacuum and it can 
be given the following definition 

2(k2) _ (g(Q)v 2 g(Q)) _ A 2 (9) 

- (q(0)q(0)) = q ■ (9) 

In our study [T7], presented here, we use the value = 0.4 GeV 2 , which is supported by several 
analyses, though values within the interval [0.35-^0.45] GeV 2 are still acceptable (sec [10,20,21] 
and references cited therein). 



1 In this work we use the "fixed-point" gauge z^A^ — 0, so that [0, z] — 1. 
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The QCD SRs with nonlocal condensates for the pion DA write |10j 

oo oo 
so 

+A$ G (x, M 2 ) + [A<f> s (x, M 2 ) + A$ v (x, M 2 ) + A$ T (x, M 2 )] q , (10) 

where (fA_ 1 is the Ai-meson DA and f n and Ja ± are, respectively, the decay constants of the 
pion and the ^i-meson. Here and below we use M 2 notation for Borel parameter. Note 
that the ^4i-meson state is an effective state which takes into account both the tt' and the 
ai meson. The nonperturbative input in the theoretical part of the SR (right side) are the 
gluon-condensate term A$>q(x, M 2 ) and the quark-condensate contribution [...]q. This latter 
contribution contains the vector-condensate term (V), the mixed quark-gluon condensate term 
(T), and the scalar condensate term (S). The explicit expressions for the perturbative and the 
nonperturbative contributions to the NLO spectral density p^llf^ (x) can be found in [TUlfTT]. 
Remarkably, the first radiative correction in the spectral density in the end-point region is of 
0{a s ) and comes out too large relative to both the zeroth-order perturbative contribution as 
well as the nonperturbative parts. For this reason, we resort in our analysis [T7j to the leading- 
order (LO) approximation {x) = 3xx/2n 2 . In fact, in order to include radiative corrections 
into the spectral density (when analyzing the end-point region), one would be obliged to resum 
all radiative corrections — a formidable task for the future. 

Among the nonperturbative terms, the scalar-quark condensate provides the largest smooth 
contribution at the origin x ~ 0, notably, 

A$ s (x,M 2 ) = 18A s {&x + <S>"x 2 + 0(x 3 )) (11) 
with the coefficients As = - ^ S (qq) 2 , 

$/ „ f°° M^l da = {qq)- 1 [°°z 2 Ms(z 2 ) dz 2 and $" « $'M 2 f s {M 2 t)— dt . (12) 
Jo « Jo Jo 1-* 

Inspection of the expression for the first coefficient $' reveals that the end-point behavior of the 
pion DA is directly related to the behavior of the scalar-quark condensate at large/moderate 
distances. 

On the other hand, all nonperturbative terms in the local condensate model (X 2 = 0) are 
concentrated exactly at the endpoints (A<]?g oc (x) ~ S(x)), as shown in Fig. [3] in terms of the 
scalar-condensate contribution. In contrast, the nonlocal condensate leads to an end-point- 
suppressed nonperturbative contribution. 

In using the SR (jlOp in order to study the end-point behavior of the pion DA, we have 
to find an appropriate characteristic that is capable of describing the slope of the pion DA at 
the origin. Because the nonperturbative terms are strongly concentrated at the endpoints (see 
Fig- E]), viz., A<3?v ~ x 5'(A - x), A$ G ~ S(A - x), . . . with A = \ 2 q /(2M 2 ) 1 the best way to 
take into account all these contributions is to investigate the integral characteristics of the pion 
DA. For this reason we invented in [T7] the "integral derivative", which will occupy us in the 
next section. 

If we apply the Gaussian model ([8]), then only the four-quark condensate A<i>s ~ x 9{A — x) 
contributes in the end-point region — see Fig. [3] — without leading to singularities. But, if 
we assume a behavior of the various condensates differing from the delta-ansatz model, and 
use, for instance, for the scalar-quark condensate a smooth model like (|2ip (which implies a 
decay at large distances not slower than the exponential decay), then the other nonperturbative 
terms (V), (G), and (T) contribute only small amounts in the small-x region. Therefore, these 
terms can be neglected if one is only interested in deriving the simplest characteristic, i.e., the 
derivative of the pion DA at the origin ^4(0), considered in Sec.[H 



5 




Figure 3: Comparison of the leading nonperturbative contribution in the local (A<J > s° c ( :z; )) an d the nonlocal 
(A$g LC (a;)) condensate models along with the perturbative contribution (p pe tt(x)). 

3 "Integral" sum rules 

The key features of the applied SRs are collected in Table Because the QCD SRs were 
developed with the aim to study the integral characteristics of the pion DA, most approaches 
appeal to the moments 

(Z N )n= f (2x-l) N p n (x)dx, (13) 
Jo 

where £ = 2x — 1. Once these moments are known, they can be used in order to reverse 
engineer the pion DA, with a precision depending upon the influence of the magnitude of the 
discarded higher-order moments. The zero-order moment leads to a SR for the decay constant, 
studied long ago by Shifman, Vainshtein and Zakharov [22J. The second-order term of the 
Gegenbauer expansion was determined by Chernyak and Zhitnitsky [12] from the standard SRs 
with local condensates (Ajj = 0). Still higher-order coefficients were computed [10] using QCD 
SRs with NLCs. It was shown there that one can de facto resort to the first two Gegenbauer 
coefficients ai and because the values of aj with i = 6,8, 10 were calculated and found to 
be negligible. On that basis, the authors of [10] obtained a bunch of two-parameter dependent 
pion DAs, among them also the BMS model, mentioned earlier, and an independent SR for the 
inverse moment ([5]). The inverse moment itself is a rather good, though not sufficient (see for 
arguments [23]), indicator for the end-point behavior of the pion DA, as we can see from the 
fifth column of Table [TJ 



Table 2: Comparison of results for the moments of the pion DA obtained within different approaches which use 
QCD SRs with local or nonlocal condensates. 



Approach 


Characteristics 


Accuracy 


Condensate 


Result 


SVZ [22] 




LO 


local 


fn 


CZ [TJ] 


ie N ),N = o, i 


LO 


local 


fir, a 2 


bms ma 


(x- 1 ), (e N ),N = 0,l,...,5 


NLO 


nonlocal 


fn, 02, «4, (z" 1 ) 


Here [H] 


¥4(0), 


LO 


nonlocal 


¥4(0), [£>Myv](*) 



In our work [T7] we probed the endpoint region of the pion DA in another way which 
employs an averaged "integral" derivative. As we explained in Sec. [21 the usefulness of these 
derivatives follows from the fact that they can be applied to QCD SRs which may even contain 
singular terms. With the help of the averaging procedure, one can take into account all these 
contributions. To expose the usefulness of this procedure, let us construct the following sequence 



6 



of average derivatives obeying the condition c/?(0) = 0: standard derivative, finite-difference 
derivative, and generalized inverse moment given by 



[D^ip](x) = ip'(x), [D^(p](x) = <p(x)/x, and [D®<p](x) = - f^-dy (14) 

x J y 



o 



respectively. The integral derivatives D^ u \ generalizing this sequence, were defined in [17] by 
means of the following expression: 



X 



[DW<p](x) = - hp(y)f(y,u-2,x)dy J(y,v,x) = -f — ^ In- at <p{0) = . (15) 

x J r(i/ + \)y v y, 



o 



When = 2 and x = 1, the integral derivative coincides with the inverse moment of the pion 
DA: [Z)( 2 )(^](l) = (x^ 1 )-^. Assuming that the Taylor expansion of f(x) at x = exists, one 
finds from (1151) 



[D^<p](x) = <p'(0) + ^(0)^|zr + O (^ (3) (0) 



(16) 



which is valid for any real v, as we explained in detail in [T7]. From the above equation, one 
can see that the defined operator reproduces at small x and/or large v the derivative of 
(p(x) at the origin x = 0. 

By applying the operator [_D("+ 2 )] on both sides of the QCD SR given by (|10p . we obtain a 
new SR for [D^+^ip^ix), viz., 

oo 

fl [£><V](*) + e- m V M2 [DM Wl ](x) + / (*) e"^^ 

(17) 



[^ H P P ert] (x)e- s / A/ + [D^A$ G ](x,M 2 ) + [L)MA$v](x,M 2 ) 

+ [L>^A$ T ](x,M 2 ) + [D^A$ s ](x,M 2 ). 

The detailed analysis of this SR can be found in [TTj. Here we only discuss the nonperturbative 
terms depicted in Fig. The dominant contribution to the integral derivative [D^ <p w ] (0.5) of 
the pion DA stems from the scalar-quark condensate (S), while the vector-condensate (V), the 
mixed quark-gluon condensate (T), and the gluon-condensate (GG) contributions are compar- 
atively less important. It turns out that the image of the operator for v > 6 is numerically 
very close to the result obtained with the differentiation method (see next section) — for any 
x. 

For large v, the (S)-term dominates and is close to the value obtained by the standard 
derivative illustrated in Fig. S]by the horizontal line, while all other condensate terms disappear. 
Thus, the integral SR (|17p becomes close to the differential one to be considered in the next 
section. For this reason, we analyze the constructed SR (|17|) for u 6 [2, 6] and x > 0.4 and 
present the results in Fig. [5] in terms of a solid line that is inside the light gray strip bounded 
by the short-dashed lines. For the sake of comparison, the predictions for the asymptotic DA 
(dashed-dotted line) and the BMS DA bunch — obtained in the NLC SR analysis of Ref. [10] 
- (shaded band limited by long-dashed lines) are also shown. From this figure we see that our 
SR estimates for [D^ (p^ R ](x) agree fairly well with the BMS model — see also Table O This 
table shows estimates for the third-order integral derivative of the pion DA for x = 0.5, using 
(i) the sum rule given by Eq. (|17p and (ii) the pion DA models we discussed above; in addition, 
results pertaining to flat-type DAs are also included. 
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Figure 4: Mutual comparison of the nonperturbative contributions (k=S, V, T, and GG) to the SR {T7J) for 
the integral derivatives [D^'ip^x) of the pion DA and with the value 18A S $'//^ = 12A s f~ 2 \~ i (horizontal 
solid line) of the scalar-quark contribution to the SR (|19[) for the standard derivative v?jt(0). The scalar-quark 
condensate term (k=S) is shown as a dashed (blue) line. The other terms are the vector-quark condensate term 
(k=V) — dashed-dotted (green) line, the mixed quark-gluon condensate term (k=T) — dashed-dotted-dotted 
(red) line, and the gluon-condensate term (k=GG) — dotted (black) line. 



First, we compare the QCD SR result, obtained from (|17p . with what one finds with the 
flat-type DA models. Consider first the flat-top DA model defined by 

= r ^["^W -*r- (18) 

This model was invented in [13J and attempts to describe the BaBar data via a logarithmic 
behavior with Q 2 . For the value a = 0.1, one finds [.£>( 3 )<^ flat ](0.5) = 227, which is much larger 
and far outside the range of values extracted from our SR flTTJ) . 

As a second option, we consider a particular flat-type pion DA which is provided by the 
AdS/QCD correspondence in the holographic approach — see, for instance, Refs. j24}426j. In 
that case one has a = 0.5 yielding [D^(p hol ](0.5) = 14. 

On the other hand, the CZ model — being also endpoint-enhanced — yields third-order in- 
tegral derivatives which are incompatible with the values derived from our SR (|17p (cf. Table [3|). 
Note that a similar statement also applies to the pion DA proposed in [27], which employs a 
Brodsky-Huang-Lepage ansatz for the k^-dependence of the pion wave function — see Table [3l 
The main message from this table is that the SR for the integral derivative of the pion DA is 



f 
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Figure 5: ^-dependence (left panel) and //-dependence (right panel) of [D^ ip v ](x) shown for the BMS bunch of 
pion DAs [TD] (shaded green band within long-dashed lines) in comparison with the SR result (|17p (narrow gray 
strip) in both panels. The left panel shows the predictions for [D^ ifi n ](x), whereas those for [_D^^ W ](0.5) are 
presented in the right panel. The dashed-dotted line denotes the asymptotic result [D^"' tp Asy ](x) — 6 — 3x/2"~ 2 . 
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fulfilled by the BMS bunch, whereas flat-type DAs and the CZ model DA have no overlap with 
the estimated range of values. By contrast, the pion DA model proposed in [27] — though it 
provides a similarly large integral derivative like the CZ DA — has a usual derivative at the 
origin which is zero due to the strong exponential suppression of this DA in the small vicinity 
of the origin. Finally, the model DA defined by Eq. (|18j) has no derivative at the origin while 
the integral derivative is well-defined. As we see from these examples, the integral derivative 
allows one to compare a broader range of the pion DA values in the end-point region than the 
standard one. 



Table 3: Results for the third-order integral derivative at the value x = 0.5 and such for the standard derivative 
of the pion DA, using different SR approaches (first three rows) and pion DA models (six last rows). 





Approach /Model 


"Integral derivative" [D^Lp n ](0.5) 


Derivative ip' n (0) 


1 


Integral SR ([17) 


4.7 ±0.5 


5.5 ± 1.5 


2 


Differential SR flU 




5.3 ±0.5 


3 


SR (USD with smooth NLC ([21]) 




7.0 ±0.7 


4 


BMS bunch [TU] 


5.7 ± 1.0 


1.7 ±5.3 


5 


Asymptotic DA 


5.25 


6 


6 


CZ DA [28] 


15.1 


26.2 


7 


DA from [27] 


14 





8 


AdS/QCD DA [2] 


14 


» 6 


9 


flat-top DA (Eq. ([15]), a = 0.1) 


227 


> 6 



It is worth mentioning that the usual derivative <p' n (0) of the pion DA encapsulates the key 
characteristics of the pion DA at small x. This quantity can be extracted from [D^(p^](x), 
which can be derived from the SR (j!7|) employing different values of v and x. The results 
are shown in Fig. [5] To determine tp' w (0) one can first use two terms of the Taylor expansion 
([15]) and then subtract the second derivative for which the asymptotic value tp"(Q) = —12(6) is 
taken. That yields the estimate <p' n (0) = 5.5 ± 1.5 for any 0.4 < x and 2 < v < 6. 

4 Differential sum rules 

Another way to study the behavior of the pion DA in the small- x region is provided by the 
differentiation of the SR (|10p . which yields 

/^(0,M 2 ) = JL M 2 _ e-,/* 9 ) + 18 As& - f Al ^ x (0)e-"V^ 2 . (19) 

We shall evaluate this SR for the threshold value sq = 2.61 GeV 2 , recalling that we are employing 
a LO expression for the spectral density. As it was shown in our recent work [IT], only the 
four-quark condensate survives and gives a contribution to the SR defined by Eq. (fl~2|) . ^From 
this equation we see that the nonperturbative contribution to the SR is mainly due to the 
scalar-quark condensate at large/moderate distances z 2 > 4/ (A; 2 ). Employing the delta-ansatz 
fs((x) = 5(a — A 2 /2) the nonperturbative contribution to SR (fl~9|) reduces to the following simple 
expression 

= delta = 4- (20) 
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Having fixed the ingredients of the SR, we are now able to consider the implications of 
using the smooth model for the quark- virtuality distribution in the differential SR. Despite the 
usefulness of the Gaussian model, we have to take into account the possibility that the scalar- 
quark condensate may behave differently at asymptotically large distances. Indeed, there are 
indications from heavy-quark effective theory [29] that it could decay exponentially. Note that 
in order to ensure the existence of the vacuum matrix element (q(D 2 ) N q) , the quark-virtuality 
distribution fs(a) should decrease faster than any power l/a N+l as a — > oo [15]. Following 
this reasoning, a two-tier model for fs was proposed in [21,30j which has a smooth dependence 
on the quark virtuality a. Hence, one has 

where K n {z) is the modified Bessel function. This so-called "smooth model", depends on two 
parameters A and a which serve to take into account the long- and the short-distance behavior 
of the nonlocal condensates. For large distances | z\ = y/—z 2 , this model leads to the asymptotic 
form of the scalar quark NLC 

Ms(z>) "±5° m^We-W ^T^ ■ (22) 

VAK n (2Aa) 

Let us now consider briefly this model — referring for technical details to [21|,I30| — and 
set n = 1, whereas the second parameter A = 0.45 GeV can be taken from the QCD SRs for 
the heavy-light meson transition in heavy quark effective theory [291 131] . The two parameters 
n and A are responsible for the large- z behavior of the scalar-quark condensate, cf. Eq. (|22|) . 



The third parameter a 2 = 10 GeV 2 is defined in terms of the parameters n, A, and A 2 via the 
following equation 



f 

Jo 



°° t ( a u A ^ n+1 (2Aa) A 2 
af s ( a ;A,n,a)da = - = - , (23) 



which we are going to evaluate for the value of the nonlocality parameter X 2 = 0.4 GeV 2 . The 
main effect of using a smooth model for the quark-virtuality distribution in comparison to the 
Gaussian form, fs(a) = 5 (a — A 2 /2), is the increase of the nonperturbative contribution to the 
SR, induced this way, entailing the relation 

/•°° / g (q;A,n,(T) a 2 K n _ 2 (2ko) 

^ooth = J Q ~ 2 da = F Kn{2A(j) > $ delta . (24) 



We studied the SR (|T9|) for this model adopting the following values of its parameters: 
f s (a; A = 0.45 GeV, n = 1, a 2 = 10 GeV" 2 ). It turns out that the average value of the derivative 
ip' n (0,M 2 ) in the fiducial Borel interval is ^(0) = 7.0(7), meaning that the nonperturbative 
contribution $g th , obtained from the smooth model, is approximately two times larger than 
the analogous contribution ^^eita f° r ^ ne delta ansatz $g mooth ~ 2.3 3>^ clta . Appealing to the 
relation (|12|) . it seems reasonable to conclude that choosing a model for the condensate that has 
a slower decay at large distances (small n or A), may induce an increase of the nonperturbative 
contribution to the SR (|10p and, hence, entail an increase of the value <^4(0) as well. The option 
of having a condensate model with a faster decay at large distances (large n or A), leads to a 
decrease of the nonperturbative contribution to the SR (fTU|) and therefore to a decrease of the 
value <^4(0). To facilitate the comparison of these two distinct possibilities for the scalar-quark 
condensate, we give in the last column of Table [3] the values of the (usual) pion DA derivative 
at x ~ 0, using various SR approaches (first three rows) and selected pion DA models (last six 
rows). 
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5 Conclusions 



In the present work we proposed a direct way to access the end-point characteristics of the pion 
DA in the QCD SR approach with nonlocal condensates. To characterize the slope of the pion 
DA at the origin, we introduced [T7] a suitable operator for integral derivatives that allows us to 
describe a broader range of the pion DA in the end-point region than the standard derivative. 
Moreover this operator provides the possibility to include all condensate terms of the QCD SRs 
with nonlocal condensates. 

Our results have been presented in the second and the third column of Table for the range 
of values of the integral and the standard derivatives of the pion DA, respectively. The same 
table contains also the values of the derivatives of some characteristic pion DA models, viz., 
the BMS, the CZ, the asymptotic, and the flat-top DA given by Eq. (JIBJ) with a = 0.1. The 
dependence of the integral derivative on its parameters x and v is displayed in Fig. [5j 

From Table El we conclude that the differential (fT9|) and the integral (fTTj) SRs agree rather 
well with each other. The integral and the standard derivatives of the pion DA, based on these 
new SRs, give smaller values than the asymptotic one and overlapping with the range of values 
determined with the BMS bunch of pion DAs |10| . while there is no agreement with the CZ DA 
and the considered flat-top model. 

It is worth remarking that, employing the integral and the differential QCD SR (|10p — which 
employ the minimal Gaussian model for the nonlocal condensates — cannot be satisfied by flat- 
type pion distribution amplitudes. Using a physically motivated exponential decay model by 
means of expression (j2Tj) . leads to a higher value of the slope of the pion DA at the origin ip' n (0), 
though it is still much smaller than the corresponding value of the flat-type pion distribution 
amplitudes. 

Some final comments: It turns out that the nonperturbative content in the differential 
SR for </4(0) is mainly due to the scalar-quark condensate, a feature valid also for the slope 
obtained with the integral SR. To be specific, the scalar-quark condensate term is proportional 
to the second inverse moment (|12p of the distribution fsip) of the vacuum-quark virtuality 
and is determined by the behavior of the quark condensate at large/moderate distances of the 
vacuum quarks. By virtue of (|12p . we may conclude that, adopting a model for the scalar-quark 
condensate that has a slower decay at large distances, entails an increase of the nonperturbative 
contribution to the SRs (fTTl and fT9|) . so that also the value of the pion DA slope defined via the 
integral and the standard derivatives increases. 

The presented analysis shows that it is difficult to reconcile QCD sum rules, and related 
techniques and features of the QCD nonperturbative vacuum in terms of nonlocal condensates, 
with flat-type pion DAs. 
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